A Lattice-Spin Mechanism in Colossal Magnetoresistant Manganites 
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We present a single-orbital double-exchange model, coupled with cooperative phonons (the so 
called breathing- modes of the oxygen octahedra in manganites). The model is studied with Monte 
Carlo simulations. For a finite range of doping and coupling constants, a first-order Metal-Insulator 
phase transition is found, that coincides with the Paramagnetic-Ferromagnetic phase transition. The 
insulating state is due to the self-trapping of every carrier within an oxygen octahedron distortion. 
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Mixed-valence manganites have attracted much atten- 
tion lately, since they undergo a transition from a ferro- 
magnetic (F) to a paramagnetic (P) state accompanied 
by a metal (M) to insulator (I) transition. The Double 
Exchange (DE) mechanism H plays a major role to ex- 
plain the magnetic transition, but the mechanism respon- 
sible of the M-I transition is not fully understood S . In 
the DE model, due to a very strong Hund's coupling, the 
carriers are strongly ferromagnetically coupled to the Mn 
core spins producing a modulation of the hopping ampli- 
tude between Mn ions. Recently a growing number of 
experiments [^-^| support the idea that phase separation 
is important in manganites |f] ^) . In this state, a phase is 
metallic and the other insulating. The theoretical chal- 
lenge is thus to find a FM — PI phase transition of the 
first order. Whereas it is widely assumed that the metal- 
lic phase is the ferromagnetic DE phase, it is not clear 
what is the origin of the gap in the insulating phase. For 
doping level near half filling and low critical temperature 
(T) materials, it has been proposed and widely accepted 
that the insulator is a charge-orbital ordered phase with 
a strong short range order. At lower doping levels the in- 
sulating phase is supposed to be some kind of polaron gas 
0,0 (polaron meaning a lattice-spin object) or polaron 
lattice phase |p^ , |l3| . 

An attractive picture for polaron formation was pre- 
sented in Refs. Jl^,|ll[: it was observed that the impor- 
tance of the electron-lattice coupling is given by its ratio 
with the carriers kinetic energy. Since within the DE 
mechanism the kinetic energy decreases upon heating, it 
was proposed that localized polarons are formed at the 
F-P transition giving rise to a M-I transition. Millis et al. 
JiT| used the dynamical mean field method to study the 
coupling of the carriers to local Jahn- Teller distortions 
and to the Mn core spins. At half filling an I-M transi- 
tion was found close to the Curie T . It was shown that 
the electron-phonon coupling could be tuned to repro- 
duce the T-dependence of the resistivity of several man- 
ganites. Unfortunately, this approach presented several 
caveats: the M-I transition was only found at half-filling, 
phonons are treated classically, and (most important) in- 



tersite phonon correlations were not considered. 

In this Letter, we consider a single-orbital (s-wave) DE 
model coupled with phonons. The model will be kept as 
simple as possible, since our scope is to shed some light 
on the mechanism behind the coupling between the M-I 
transition and the Curie temperature. Both core-spins 
and phonons arc treated as classical variables (for spins 
this is a controlled approximation H]). The lattice dis- 
tortion we study is the deformation of the oxygen octahe- 
dra around Mn sites. The coupling of these modes with 
charge carriers is expected to be at least as large as the 
one producing Jahn- Teller distortions [Gl| . The contrac- 
tion {breathing) of a MnOg octahedron, implies a volume 
growth of its neighbors (but do not change the total lat- 
tice volume, as would be needed to study magnetostric- 
tion effects |l6[|). Thus this mode is strongly cooperative, 
and not suitable for Mean-Field studies. A Monte Carlo 
(MC) investigation is therefore performed. The super- 
exchange antiferromagnetic coupling between the core 
spins [0 is neglected in our model. Also, as we men- 
tioned above, we work with a single s-orbital per site and 
we do not consider the two degenerated e g orbitals which 
are crucial to understand the magnetic phase-diagram 
beyond half- filling |H. Thus, our model is to be re- 
garded as a model for materials like Lai-^Ca^MnOa in 
the 0.15 < x < 0.4 regime where the magnetoresistance 
is largest, and the only experimentally relevant magnetic 
phases are F and P |. In spite of its simplicity, the 
model presents a first-order M-I phase transition, that 
coincides with the P-F phase transition for a finite-range 
of doping and coupling constants, in contrast with previ- 
ous work ]l0| . 

The model Hamiltonian contains the Mn e g itinerant 
carriers coupled to the t g Mn core spins, and phonons: 

H = i?KE + #Hund + -^cl-ph + -ffph ■ (1) 

Here Hke is the kinetic energy of the carriers hopping be- 
tween Mn atoms that form a simple-cubic lattice, -fffjund 
is the Hund interaction, i/ i_ p h is the lattice-carrier cou- 
pling energy, while -ff p h represents the crystal elastic 
energy. The Hund interaction is very large in man- 
ganese oxides and at each place the carrier spin is forced 
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to be parallel to the core spin, which allows to reduce 
Hke + ^Hund to the DE Hamiltonian M : 
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T{Si,Sj) c\cj + h 



(2) 



Here cf creates an electron at place i with spin parallel 
to the core spin at z, t is the hopping amplitude between 
first neighbors ions and T(Si,Sj) = — t [cos ^ cos + 
sin 4f sin %^ >i ^' p ^\ 9i and (pi being the polar coordi- 
nates of the core spin at site i, Si. For the electron- 
phonon coupling we consider the distortions of the MnC>6 
octahedron formed by the six oxygens surrounding the 
Mn ions. The oxygens are located at the center of the 
edges of the cubic lattice formed by Mn atoms. Each 
oxygen is allowed to move along the edge on which it 
is located. The distortions of the six oxygens surround- 
ing a Mn at site i are given by Ui t ± a where a run over 
x, y, z. The size fluctuations of the MnC-6 octahedra 
are coupled to charge fluctuations in the Mn through the 
electron-phonon interaction, 



(3) 



where A is the electron-phonon coupling. This interaction 
tends to produce lattice distortions. This tendency is 
opposed by the stiffness of the Mn-0 bonds: 
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FIG. 1. The T (in units of i) — A phase diagram of the 
model ([!]), obtained from MC simulations on a 6 3 lattice with 
x — 0.08 (see text). Filled dots: P-F critical T. Open trian- 
gles: M-I critical T. 



One can get some intuition about the physics of our 
model considering the limit of very few carriers. If the 
carrier-lattice coupling is strong, one can gain enough 
electronic energy by contracting an oxygen octahedron 
(thus localizing a carrier: a polaron) to compensate the 
high price in elastic energy. For the fully spin-polarized 
lattice, one finds (using e.g. the techniques of Ref. |19|) 



^threshold _ ^ 91 The carrier is localized at the polaron's 
center with a 96% probability. If one repeats the calcu- 
lation for the P phase, using deGennes' virtual-crystal 
approximation @, finds A tllrcshold » 1.56 (the P po- 
laron is also at its center with 96% probability). Thus 
for A in between both thresholds, one expects that at the 
F-P phase transition, every carrier will form a strongly 
localized polaron upon heating. The system becomes in- 
sulating due to the formation of a fully occupied band 
separated from upper states by a gap. This picture is 
largely confirmed by the MC simulation of the model. 
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FIG. 2. Top: M vs. T for A = 0, 1.35, 1.45 and 1.7, for 
the MC simulations in Fig. |l| Bottom: (decimal) logarithm 
of the average d.c. resistance (a.u.) vs. T, from the same 
simulations (A = 1.7 data are out range). 

For calculating the T-dependent phase diagrams, we 
perform MC simulations on the classical variables: the 
core spins Si and the oxygen displacements {ui^ a ). The 
simulations are done in N x N x N lattices with periodic 
boundary conditions, using an standard Metropolis algo- 
rithm. The kinetic energy of the carriers is calculated by 
diagonalizing the electron Hamiltonian at each Metropo- 
lis step. The diagonalization CPU cost grows like N 6 
and has limited us to N — 6 pl| . We have used the 
N — 4 results to check for finite-size effects. Altghough 
both absolute values and details sometimes change for 
N — 4 simulations, all issues discussed in this letter re- 
main valid. The Fermi temperature of the carriers is 
much higher than other T in the system and we assume 
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the carriers to be at zero T p3| . We calculate the ther- 
mal average of different physical quantities; the absolute 
value of the Si polarization, M, the electronic energy dif- 
ference between the lowest energy empty state and the 
highest energy occupied one, E gap , the standard devia- 
tion of the (spatial) probability distribution of the MnOg 
octahedra volume, AV rms , and the electronic density of 
states, p(uj). We also measure the average d.c. resistance 
of the system, calculating the resistance of the N x N x N 
cubic lattice connected to two semi-infinite perfect leads 
p4| , p5[ using the standard Kubo formula pc| . 
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FIG. 3. Standard deviation of tne (lattice) distribution of 
the volume of the MnOe octahedra vs. T, for the simulations 
in Fig. | 

Fig. ^ shows the phase diagram, A vs. T, for a N =6 
cubic lattice with 17 electrons, i.e. x ~0.08. We use the 
criterion that the system is M/I when the d.c. resistance 
increases/decreases with T. The phase diagram contains 
four phases: FM, FI, PM and PI. As expected for small A, 
polarons are not formed by the (small) lattice distortions, 
the system being metallic at all T [^3],^5| (see the growing 
behavior of the d.c. resistance upon heating in Fig. |^ — 
bottom). When T grows there is a second order FM-PM 
transition, (see the smooth temperature behavior of M 
in the top of Fig. ||). At A=0, we found the Curie T 
at 0.0722;, in agreement with previous MC simulations at 
£=0.08 ^3|. Notice that thermodynamic quantities are 
even functions of A (because of the symmetry A — » —A, 
Ui.a —> —Ui,a hi E<l-(0))i and thus for small couplings 
they quadratically depend on A. This can be checked for 
the Curie T in Fig. |. 



At intermediate A, new behavior is expected. In the 
P phase (that has higher energy than the F phase), all 
carriers form polarons and the system is an insulator, 
while in the F phase there are not polarons and the sys- 
tem is metallic. Upon heating, the degenerated electronic 
system undergoes a phase transition at the F-P transi- 
tion, with a sharp change in electronic energy. A first- 
order F-P phase transition with growing T is obtained, 
that coincides with a M-I transition. Indeed, Fig. [j] 
shows how the P-M and I-M transition lines merge for 
1.4 < A < 1.65. M changes abruptly at the phase tran- 
sition (Fig. ||, top) while the d.c. resistance grows by 
a factor 10 6 (Fig. |[ bottom) and becomes a decreasing 
function of T. First-order P-F transitions are experimen- 
tally found in manganites as the strength of the electron- 
phonon coupling increases [p7| . For slightly smaller val- 
ues of A (A = 1.35 < A c = 1.4) the d.c. resistance al- 
ways grows with T, the system being M, and the F-P 
phase transition is continuous (Fig. ^, top and bottom). 
Polaron formation can also be seen in the spatial distri- 
bution of volumes of the MnOg octahedra that becomes 
very inhomogeneous: on polarons, octahedra are small 
while in most sites the volume is uniform. Note the sharp 
change of Al^ms at the critical line for 1.4 < A < 1.7 
(Fig. ||), and the smoother T-dependence for A = 1.35. 
Notice also (Fig. |) the gap in p(u) for 1.4 < A < 1.7 
(the Fermi level is in the gap). 
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FIG. 4. Thermal average of the electronic density of states 
for A = 1.45, T = 0.061 (PI phase on Fig. |). The dashed 
line is the density of states of the fully spin-polarized system, 
without octahedra distortions. 

For large A, polarons exist in both the P and F phases. 
The F-P and the M-I transitions decouple for A > 1.65 
(Fig. |l|). The F-P one is again continuous: M evolves 
smoothly with T (Fig. [| top), and the system has po- 
larons at the lowest T (Fig. ||), being an I. The F phase 
in a DE system with a fully occupied band is somehow 
unconventional, because the mean carriers energy is at 
the band center which is usually spin independent. This 
is not the case for our model. The difference between the 
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positions of the (polaronic) band center of the fully polar- 
ized and unpolarized systems can be calculated as before, 
with deGennes' virtual-crystal approximation. The en- 
ergy difference for large A is roughly — 1.08t/A 2 which 
explains the (small) ferromagnetic interaction that de- 
creases with increasing A. 

The results shown up to now can be understood qual- 
itatively and semiquantitatively within the few carriers 
limit. However, when the polaron density approaches the 
percolation threshold of the cubic lattice (p c = 0.31 ps|), 
the polaronic wave function becomes far less localized. In 
Fig. ^[ we show the phase diagram for x — 0.3. We find 
the same phases as in the x = 0.08 case. In the I case, all 
carriers are polarons. Notice the absence of a coupling 
between the P-M and the M-I transitions, the former be- 
ing of the second order. Therefore in our model, we do 
not get a T-dependent M-I transition at x — 0.3. Nev- 
ertheless our model only includes an orbital per Mn ion. 
The use of a more realistic band structure, will increase 
the phase space for the polarons and a T-dependent M-I 
transition will occur also at higher doping levels. 
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FIG. 5. Phase diagram for the n&del (§), as obtained from 
the MC simulation on a 4 3 lattices at x = 0.29. Filled dots: 
the F-P critical T. Open-triangles: M-I critical T. 

In summary, we have studied a DE model coupled with 
the breathing modes of the MnOg octahedra. Given the 
collective nature of these modes, the spatial distribution 
of the lattice distortions is inhomogeneous and a Mean 
Field study difficult. A throughout MC investigation of 
the phase diagram has been carried out. Four phases 
have been found: FM, FI, PM and PI. A first order M-I 
transition with growing T (coincident with the F-M tran- 
sition) has been found for the first time on a MC simula- 
tion. This transition survives for a finite range of doping 
and coupling constant. The M-I transition is induced by 
the self-trapping of every carrier on a polaron, making the 
Fermi level to lie on the gap between the polaronic and 
DE bands (for t » 0.16eV Bj, the gap will be optical). 
We argue that the mechanism presented here is relevant 
for the formation of the paramagnetic-insulating phase 



experimentally observed in the phase-separated state of 
colossal magnetoresistive manganites || || . 
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